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Detrended fluctuation analysis (DFA) is a scaling analysis method used to quantify long-range 
power-law correlations in signals. Many physical and biological signals are "noisy", heterogeneous 
and exhibit different types of nonstationarities, which can affect the correlation properties of these 
signals. We systematically study the effects of three types of nonstationarities often encountered 
in real data. Specifically, we consider nonstationary sequences formed in three ways: (i) stitching 
together segments of data obtained from discontinuous experimental recordings, or removing some 
noisy and unreliable parts from continuous recordings and stitching together the remaining parts 
— a "cutting" procedure commonly used in preparing data prior to signal analysis; (ii) adding to a 
signal with known correlations a tunable concentration of random outliers or spikes with different 
amplitude, and (iii) generating a signal comprised of segments with different properties — e.g. 
different standard deviations or different correlation exponents. We compare the difference between 
the scaling results obtained for stationary correlated signals and correlated signals with these three 
types of nonstationarities. We find that introducing nonstationarities to stationary correlated signals 
leads to the appearance of crossovers in the scaling behavior and we study how the characteristics 
of these crossovers depend on: (a) the fraction and size of the parts cut out from the signal; (b) the 
concentration of spikes and their amplitudes; (c) the proportion between segments with different 
standard deviations or different correlations; and (d) the correlation properties of the stationary 
signal. We show how to develop strategies for pre-processing "raw" data prior to analysis, which 
will minimize the effects of nonstationarities on the scaling properties of the data and how to interpret 
the results of DFA for complex signals with different local characteristics. 



I. INTRODUCTION 



In recent years, there has been growing evidence indi- 
cating that many physical and biological systems have no 
characteristic length scale and exhibit long-range power- 
law correlations. Traditional approaches such as the 
power-spectrum and correlation analysis are suited to 
quantify correlations in stationary signals . However, 
many signals which are outputs of complex physical and 
biological systems are nonstationary — the mean, stan- 
dard deviation and higher moments, or the correlation 
functions are not invariant under time translation 0,^. 
Nonstationarity, an important aspect of complex vari- 
ability, can often be associated with different trends in 
the signal or heterogeneous segments (patches) with dif- 
ferent local statistical properties. To address this prob- 
lem, detrended fluctuation analysis (DFA) was developed 
to accurately quantify long-range power-law correlations 
embedded in a nonstationary time series Q . This 
method provides a single quantitative parameter — the 
scaling exponent a — to quantify the correlation prop- 
erties of a signal. One advantage of the DFA method 
is that it allows the detection of long-range power-law 
correlations in noisy signals with embedded polynomial 
trends that can mask the true correlations in the fluctua- 
tions of a signal. The DFA method has been successfully 
applied to research fields such as DNA p|-p^ , cardiac 
dynamics p7|-^, human gait [38|, meteorology cli- 
mate temperature fluctuations 1 4C , river flow and dis- 
charge neural receptors in biological systems p5| . 



and economics |4q-|58[ . The DFA method may also help 
identify different states of the same system with dif- 
ferent scaling behavior — e.g., the scaling exponent a 
for heart-beat intervals is different for healthy and sick 
individuals] 1^, ^ as well as for waking and sleeping 
states ||,|3§1. 

To understand the intrinsic dynamics of a given sys- 
tem, it is important to analyze and correctly interpret 
its output signals. One of the common challenges is that 
the scaling exponent is not always constant (indepen- 
dent of scale) and crossovers often exist — i.e., the value 
of the scaling exponent a differs for different ranges of 
scales |^|5^, |o|. A crossover is usually due to a 

change in the correlation properties of the signal at dif- 
ferent time or space scales, though it can also be a result 
of nonstationarities in the signal. A recent work consid- 
ered different types of nonstationarities associated with 
different trends (e.g., polynomial, sinusoidal and power- 
law trends) and systematically studied their effect on 
the scaling behavior of long-range correlated signals . 
Here we consider the effects of three other types of non- 
stationarities which are often encountered in real data or 
result from "standard" data pre-processing approaches. 

(i) Signals with segments removed 

First we consider a type of nonstationarity caused by dis- 
continuities in signals. Discontinuities may arise from the 
nature of experimental recordings - e.g., stock exchange 
data are not recorded during the nights, weekends and 
holidays p6|~p3[. Alternatively, discontinuities may be 
caused by the fact that some noisy and unreliable por- 
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tions of continuous recordings must be discarded, as often 
occurs when analyzing physiological signals ||l7|-p7|. In 
this case, a common pre-processing procedure is to cut 
out the noisy, unreliable parts of the recording and stitch 
together the remaining informative segments before any 
statistical analysis is performed. One immediate prob- 
lem is how such cutting procedure will affect the scaling 
properties of long-range correlated signals. A careful 
consideration should be given when interpreting results 
obtained from scaling analysis, so that an accurate es- 
timate of the true correlation properties of the original 
signal may be obtained. 

(ii) Signals with random spikes 

A second type of nonstationarity is due to the existence 
of spikes in data, which is very common in real life sig- 
nals jl^-Q. Spikes may arise from external conditions 
which have little to do with the intrinsic dynamics of 
the system. In this case, we must distinguish the spikes 
from normal intrinsic fluctuations in the system's output 
and filter them out when we attempt to quantify corre- 
lations. Alternatively, spikes may arise from the intrinsic 
dynamics of the system, rather than being an epiphe- 
nomenon of external conditions. In this second case, 
careful considerations should be given as to whether the 
spikes should be filtered out when estimating correla- 
tions in the signal, since such "intrinsic" spikes may be 
related to the properties of the noisy fluctuations. Here, 
we consider only the simpler case - namely, when the 
spikes are independent of the fluctuations in the signal. 
The problem is how spikes affect the scaling behavior 
of correlated signals, e.g., what kind of crossovers they 
may possibly cause. We also demonstrate to what ex- 
tent features of the crossovers depend on the statistical 
properties of the spikes. Furthermore, we show how to 
recognize if a crossover indeed indicates a transition from 
one type of underlying correlations to a different type, 
or if the crossover is due to spikes without any transition 
in the dynamical properties of the fluctuations. 

{Hi) Signals with different local behavior 
A third type of nonstationarity is associated with the 
presence of segments in a signal which exhibit different 
local statistical properties, e.g., different local standard 
deviations or different local correlations. Some examples 
include: (a) 24 hour records of heart rate fluctuations are 
characterized by segments with larger standard deviation 
during stress and physical activity and segments with 
smaller standard deviation during rest |l9[ ; (b) studies 
of DNA show that coding and non-coding regions are 
characterized by different types of correlations [|[ || ; (c) 
brain wave analysis of different sleep stages (rapid eye 
movement [REM] sleep, light sleep and deep sleep) indi- 
cates that the signal during each stage may have different 
correlation properties p^ ; (d) heartbeat signals during 
different sleep stages exhibit different scaling proper- 
ties Q. For such complex signals, results from scaling 
analysis often reveal a very complicated structure. It is a 



challenge to quantify the correlation properties of these 
signals. Here, we take a first step toward understanding 
the scaling behavior of such signals. 

We study these three types of nonstationarities embed- 
ded in correlated signals. We apply the DFA method to 
stationary correlated signals and identical signals with 
artificially imposed nonstationarities, and compare the 
difference in the scaling results, (i) We find that cutting 
segments from a signal and stitching together the remain- 
ing parts does not affect the scaling for positively corre- 
lated signals. However, this cutting procedure strongly 
affects anti-correlated signals, leading to a crossover from 
an anti-correlated regime (at small scales) to an uncor- 
related regime (at large scales), (ii) For the correlated 
signals with superposed random spikes, we find that the 
scaling behavior is a superposition of the scaling of the 
signal and the apparent scaling of the spikes. We analyt- 
ically prove this superposition relation by introducing a 
superposition rule. (Hi) For the case of complex signals 
comprised of segments with different local properties, we 
find that their scaling behavior is a superposition of the 
scaling of the different components — each component 
containing only the segments exhibiting identical statis- 
tical properties. Thus, to obtain the scaling properties 
of the signal, we need only to examine the properties of 
each component — a much simpler task than analyzing 
the original signal. 

The layout of the paper is as follows: In Sec. ||, we 
describe how we generate signals with desired long-range 
correlation properties and introduce the DFA method to 
quantify these correlations. In Sec. IH, we compare the 



scaling properties of correlated signals before and after 
removing some segments from the signals. In Sec. EV, we 



consider the effect of random spikes on correlated signals. 
We show that the superposition of spikes and signals can 
be explained by a superposition rule derived in Appendix 
In Sec. 0, we study signals comprised of segments 
with different local behavior. We systematically examine 
all resulting crossovers, their conditions of existence, and 
their typical characteristics associated with the different 
types of nonstationarity. We summarize our findings in 
Sec. fvl. 



II. METHOD 

Using a modified Fourier filtering method | |63[ |, we 
generate stationary uncorrelated, correlated, and anti- 
correlated signals u{i) (i = 1,2,3,..., A^max) with a stan- 
dard deviation cr = 1. This method consists of the fol- 
lowing steps: 

(a) First, we generate an uncorrelated and Gaussian 
distributed sequence ri{i) and calculate the Fourier trans- 
form coefficients rj^q). 

(b) The desired signal u{i) must exhibit correlations, 
which are defined by the form of the power spectrum 



S{q) = {u{q)u{-q)) - q 



-(1-7) 



(1) 
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where u[q) are the Fourier transform coefRcients of u{i) 
and 7 is the correlation exponent. Thus, we generate 
u{q) using the following transformation: 



uiq) ^ [5(g)]l/^((z), 



(2) 



where S{q) is the desired power spectrum in Eq. (^. 

(c) We calculate the inverse Fourier transform of u{q) 
to obtain u{i). 

We use the stationary correlated signal u{i) to generate 
signals with different types of nonstationarity and apply 
the DFA method to quantify correlations in these non- 
stationary signals. 

Next, we briefly introduce the DFA method, which in- 
volves the following steps 

(i) Starting with a correlated signal u{i), where i = 
1, ..,Nmax and N^ax is the length of the signal, we first 
integrate the signal u{i) and obtain y{k) = ~ 
(m)], where (u) is the mean. 

(ii) The integrated signal y{k) is divided into boxes of 
equal length n. 

(in) In each box of length n, we fit j/(fc), using a poly- 
nomial function of order £ which represents the trend in 
that box. The y coordinate of the fit line in each box 
is denoted by (see Fig. |], where linear fit is used). 

Since we use a polynomial fit of order we denote the 
algorithm as 'DFA-l. 
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FIG. 1. (a) The correlated signal u{i). (b) The integrated 
signal: y{k) — ~ W)]- The vertical dotted lines in- 

dicate a box of size n — 100, the solid straight lines segments 
are the estimated linear "trend" in each box by least-squares 
fit. 



(iv) The integrated signal y{k) is detrended by sub- 
tracting the local trend yn{k) in each box of length n. 

(v) For a given box size n, the root mean-square 
(r.m.s.) fluctuation for this integrated and detrended sig- 
nal is calculated: 



F{n) 



\ 



N„ 



Nma, 

E 

fe=i 



[v{k)-y^{kW 



(3) 



(wi) The above computation is repeated for a broad range 
of scales (box sizes n) to provide a relationship between 
F{n) and the box size n. 

A power-law relation between the average root-mean- 
square fluctuation function F(n) and the box size n in- 
dicates the presence of scaling: F{n) ^ n". The fluctu- 
ations can be characterized by a scaling exponent a, a 
self-similarity parameter which represents the long-range 
power-law correlation properties of the signal. If a = 0.5, 
there is no correlation and the signal is uncorrelated 
(white noise); if a < 0.5, the signal is anti-correlated; 
if a > 0.5, the signal is correlated |64|. 

We note that for anti-correlated signals, the scaling 
exponent obtained from the DFA method overestimates 
the true correlations at small scales To avoid this 

problem, one needs first to integrate the original anti- 
correlated signal and then apply the DFA method [ |6l| . 
The correct scaling exponent can thus be obtained from 
the relation between n and F{n)/n [instead of F{n)]. In 
the following sections, we first integrate the signals under 
consideration, then apply DFA-2 to remove linear trends 
in these integrated signals. In order to provide a more 
accurate estimate of F{n), the largest box size n we use 
is iVma£c/10, where Nmax is the total number of points in 
the signal. 

We compare the results of the DFA method obtained 
from the nonstationary signals with those obtained from 
the stationary signal u(i) and examine how the scal- 
ing properties of a detrended fluctuation hmction F(n) 
change when introducing different types of nonstationar- 
ities. 



III. SIGNALS WITH SEGMENTS REMOVED 



In this section, we study the effect of nonstationar- 
ity caused by removing segments of a given length from 
a signal and stitching together the remaining parts — 
a "cutting" procedure often used in pre-processing data 
prior to analysis. To address this question, we first gen- 
erate a stationary correlated signal u{i) (see Sec. 0) of 
length Nmax and a scaling exponent a, using the mod- 
ified Fourier filtering method |6^. Next, we divide this 
signal into Nmax/W non-overlapping segments of size W 
and randomly remove some of these segments. Finally, 
we stitch together the remaining segments in the signal 
u(i) [Fig. |^(a)], thus obtaining a surrogate nonstationary 
signal which is characterized by three parameters: the 
scaling exponent a, the segment size W and the fraction 
of the signal u{i), which is removed. 
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(e) Correlated signal: a>0.5, W=20 
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FIG. 2. Effects of the "cutting" procedure on the scaling 
behavior of stationary correlated signals. Nmax = 2^" is the 
number of points in the signals (standard deviation o = 1) 
and W is the size of the cutout segments, (a) A stationary 
signal with 10% of the points removed. The removed parts 
are presented by shaded segments of size Vl^ = 20 and the 
remaining parts are stitched together, (b) Scaling behavior 
of nonstationary signals obtained from an anti-correlated sta- 
tionary signal (scaling exponent a < 0.5) after the cutting 
procedure. A crossover from anti-correlated to uncorrelated 
(a = 0.5) behavior appears at scale Ux- The crossover scale 
rix decreases with increasing the fraction of points removed 
from the signal. We determine rix based on the difference A 
between the logarithm of F{n)/n for the original stationary 
anti-correlated signal (q = 0.1) and the nonstationary signal 
with cutout segments: rix is the scale at which A > 0.04. 
Dependence of the crossover scale rix on the fraction (c) and 
on the size W (d) of the cutout segments for anti-correlated 
signals with different scaling exponent a. (e) Cutting proce- 
dure applied to correlated signals (q > 0.5). In contrast to 
(b), no discernible effect on the scaling behavior is observed 
for different values of the scaling exponent a, even when up 
to 50% of the points in the signals are removed. 

We find that the scaling behavior of such a nonstation- 
ary signal strongly depends on the scaling exponent a of 
the original stationary correlated signal u{i). As illus- 
trated in Fig. |(b), for a stationary anti-correlated signal 
with a — 0.1, the cutting procedure causes a crossover 
in the scaling behavior of the resultant nonstationary 
signal. This crossover appears even when only 1% of 
the segments are cut out. At the scales larger than the 
crossover scale rix the r.m.s. fluctuation function behaves 
as F{n) ~ n°'^, which means an uncorrelated random- 
ness, i.e., the anti-correlation has been completely de- 
stroyed in this regime. For all anti-correlated signals with 
exponent a < 0.5, we observe a similar crossover behav- 
ior. This result is surprising, since researchers often take 
for granted that a cutting procedure before analysis does 
not change the scaling properties of the original signal. 
Our simulation shows that this assumption is not true, 
at least for anti-correlated signals. 

Next, we investigate how the two parameters — the 
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segment size W and the fraction of points cut out from 
the signal — control the effect of the cutting procedure 
on the scaling behavior of anti-correlated signals. For 
the fixed size of the segments (W = 20), we find that 
the crossover scale nx decreases with increasing the frac- 
tion of the cutout segments [Fig. ||(c)]. Furthermore, for 
anti-correlated signals with small values of the scaling ex- 
ponent a, e.g., a = 0.1 and a = 0.2, we find that and 
the fraction of the cutout segments display an approxi- 
mate power-law relationship. For a fixed fraction of the 
removed segments, we find that the crossover scale 
increases with increasing the segment size W [Fig. |^(d)]. 
To minimize the effect of the cutting procedure on the 
correlation properties, it is advantageous to cut smaller 
number of segments of larger size W. Moreover, if the 
segments which need to be removed are too close (e.g., 
at a distance shorter than the size of the segments), it 
may be advantageous to cut out both the segments and 
a part of the signal between them. This will effectively 
increase the size of the segment W without substantially 
changing the fraction of the signal which is cut out, lead- 
ing to an increase in the crossover scale ■ Such strategy 
would minimize the effect of this type of nonstationarity 
on the scaling properties of data. For small values of 
the scaling exponent a (a < 0.25), we find that nx and 
W follow power-law relationships [Fig. ^(d)]. The rea- 
son we do not observe a power-law relationship between 
rix and W and between nx and the fraction of cutout 
segments for the values of the scaling exponent a close 
to 0.5 may be due to the fact that the crossover regime 
becomes broader when it separates scaling regions with 
similar exponents, thus leading to uncertainty in defining 
n X . For a fixed W and a fixed fraction of the removed 
segments [see Figs. ||(c) and (d)], we observe that Ux in- 
creases with the increasing value of the scaling exponent 
a, i.e., the effect of the cutting procedure on the scal- 
ing behavior decreases when the anti-correlations in the 
signal are weaker (a closer to 0.5). 

Finally, we consider the case of correlated signals u{i) 
with 1.5 > a > 0.5. Surprisingly, we find that the scaling 
of correlated signals is not affected by the cutting pro- 
cedure. This observation remains true independently of 
the segment size W — from very small = 5 up to very 
large W = 5000 segments — even when up to 50% of the 
segments are removed from a signal with N^ax ~ 10^ 
points [Fig. 1(e)] . 
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IV. SIGNALS WITH RANDOM SPIKES 



In this section, we consider nonstationarity related to 
the presence of random spikes in data and we study the 
effect of this type of nonstationarity on the scaling prop- 
erties of correlated signals. First, we generate surrogate 
nonstationary signals by adding random spikes to a sta- 
tionary correlated signal u(i) [see Sec. Hand Fig.|(a-c)]. 
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FIG. 3. Effects of random spikes on the scaling behav- 
ior of stationary correlated signals, (a) An example of an 
anti-correlated signal u{i) with scaling exponent a = 0.2, 
Nmax ~ 2^'^ and standard deviation a — 1. (b) A se- 
ries of uncorrelated spikes {asp = 0.5) at 5% randomly 
chosen positions (concentration p — 0.05) and with uni- 
formly distributed amplitudes Asp in the interval [—4,4]. 
(c) The superposition of the signals in (a) and (b). (d) 
Scaling behavior of an anti-correlated signal u{i) {a = 0.2) 
with spikes {Asp = 1, p = 0.05, asp — 0.5). For 
n < Hx, F{n)/n « Fr,{n)/n ~ n", where _Fl,(n)/n is 
the scaling function of the signal u{i). For n > Ux, 
F{n)/n « Fsp{n)/n ~ n""'' . (e) Scaling behavior of a cor- 
related signal u{i) {a = 0.8) with spikes {Asp = 10, p = 0.05, 
asp = 0.5). For n < nx, F{n)/n ^ Fsp{n)/n ~ 
For n > Ux, F{n)/n ^ F,,{n)/n ~ n". Note that when 
a = asp = 0.5, there is no crossover. 

We find that the correlation properties of the nonsta- 
tionary signal with spikes depend on the scaling exponent 
a of the stationary signal and the scaling exponent asp 
of the spikes. When uncorrelated spikes {ugp = 0.5) are 
added to a correlated or anti-correlated stationary signal 
[Fig ^(d) and (e)] , we observe a change in the scaling be- 
havior with a crossover at a characteristic scale riy. For 
anti-correlated signals {a < 0.5) with random spikes, we 
find that at scales smaller than Ux, the scaling behav- 
ior is close to the one observed for the stationary anti- 
correlated signal without spikes, while for scales larger 
than Hx , there is a crossover to random behavior. In the 
case of correlated signals (a > 0.5) with random spikes, 
we find a different crossover from uncorrelated behavior 
at small scales, to correlated behavior at large scales with 
an exponent close to the exponent of the original station- 
ary correlated signal. Moreover, we find that spikes with 
a very small amplitude can cause strong crossovers in 
the case of anti-correlated signals, while for correlated 
signals, identical concentrations of spikes with a much 
larger amplitude do not affect the scaling. Based on 
these findings, we conclude that uncorrelated spikes with 
a sufficiently large amplitude can affect the DFA results 
at large scales for signals with a < 0.5 and at small scales 
for signals with a > 0.5. 

To better understand the origin of this crossover be- 
havior, we first study the scaling of the spikes only [see 
Fig. ^(b)]. By varying the concentration p {0 < p < I) 
and the amplitude A^p of the spikes in the signal, we 
find that for the general case when the spikes may be 
correlated, the r.m.s. fluctuation function behaves as 

F,p{n)/n^ ko^Aspu"'", (4) 

where fco is a constant and asp is the scaling exponent of 
the spikes. 

Next, we investigate the analytical relation between 
the DFA results obtained from the original correlated sig- 
nal, the spikes and the superposition of signal and spikes. 
Since the original signal and the spikes are not correlated, 
we can use a superposition rule (see [^ and Appendix ^) 



to derive the r.m.s. fluctuation function F{n)/n for the 
correlated signal with spikes: 

[F(n)/n]2 = [F„(n)/n]2 + [Fsp{n)/n]^ (5) 

where Fr,{n)/n and Fsp{n)/n are the r.m.s. fluctua- 
tion function for the signal and the spikes, respec- 
tively. To confirm this theoretical result, we calculate 
^[F,-,{n)/nY + [Fsp{n)/nf [see Figs. |(d), (e)] and find 
this Eq. (^) is remarkably consistent with our experimen- 
tal observations. 

Using the superposition rule, we can also theoretically 
predict the crossover scale rix as the intercept between 
Frj{n)/n and Fsp{n)/n, i.e., where F^{nx) = Fsp{nx). 
We find that 

nx = ( VP^sp^j , (6) 

since the r.m.s. fiuctuation function for the signal and 
the spikes are F^(ri)/n = bon" and Fsp{n)/n — 
koy^AspH^'p [Eq. (Q)], respectively. This result predicts 
the position of the crossover depending on the parame- 
ters defining the signal and the spikes. 

Our result derived from the superposition rule can be 
useful to distinguish two cases: («) the correlated station- 
ary signal and the spikes are independent (e.g., the case 
when a correlated signal results from the intrinsic dy- 
namics of the system while the spikes are due to external 
perturbations); and {ii) the correlated stationary signal 
and the spikes are dependent (e.g., both the signal and 
the spikes arise from the intrinsic dynamics of the sys- 
tem). In the latter case, the identity in the superposition 
rule is not correct (see Appendix 



V. SIGNALS WITH DIFFERENT LOCAL 
BEHAVIOR 



Next, we study the effect of nonstationarities on com- 
plex patchy signals where different segments show differ- 
ent local behavior. This type of nonstationarity is very 
common in real world data |^,||,|l9|,|3^, Our discus- 
sion of signals composed of only two types of segments is 
limited to two simple cases: (A) different standard devi- 
ations and (B) different correlations. 
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(d) Correlated signal: a>0.5 
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FIG. 4. Scaling behavior of nonstationary correlated 
signals with diflterent local standard deviation. (a) 
Anti-correlated signal (a = 0.1) with standard deviation 
ai = 1 and amplified segments with standard deviation 
(72 = 4. The size of each segment is W = 20 and the fraction 
of the amplified segments is p = 0.1 from the total length of 
the signal {Nmax ~ 2^°). (b) Scaling behavior of the signal 
in (a) for a different fraction p of the amplified segments (af- 
ter normalization of the globe standard deviation to unity). 
A crossover from anti-correlated behavior (a = 0.1) at small 
scales to random behavior (a = 0.5) at large scales is ob- 
served, (c) Dependence of the crossover scale rix on the frac- 
tion p of amplified segments for the signal in (a). Ux is de- 
termined from the difference A of logjQ [F(n)/n] between the 
nonstationary signal with amplified segments and the origi- 
nal stationary signal. Here we choose A = 0.04. (d) Scaling 
behavior of nonstationary signals obtained from correlated 
stationary signals (1 > a > 0.5) with standard deviation 
ai — 1, for a different fraction of the amplified segments with 
a2 ~ 4. No difference in the scaling is observed, compared to 
the original stationary signal. 

Here we consider nonstationary signals comprised of 
segments with the same local scaling exponent, but dif- 
ferent local standard deviations. We first generate a sta- 
tionary correlated signal u(i) (see Sec. ||) with fixed stan- 
dard deviation ai = 1. Next, we divide the signal u{i) 
into non-overlapping segments of size W. We then ran- 
domly choose a fraction p of the segments and amplify 
the standard deviation of the signal in these segments, 
CT2 = 4 [Fig.||(a)]. Finally, we normalize the entire signal 
to global standard deviation cr = 1 by dividing the value 
of each point of the signal by ^^(1 — pjaf + pa^- 

For nonstationary anti-correlated signals {a < 0.5) 
with segments characterized by two different values of 
the standard deviation, we observe a crossover at scale 
Tlx [Fig.@(b)]. For small scales n < n^, the behavior is 
anti-correlated with an exponent equal to the scaling ex- 
ponent a of the original stationary anti-correlated signal 
u{i). For large scales n > rtx, we find a transition to 
random behavior with exponent 0.5, indicating that the 
anti-correlations have been destroyed. The dependence 
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of crossover scale rix on the fraction p of segments with 
larger standard deviation is shown in Fig. |4|(c). The de- 
pendence is not monotonic because for p = and p — I, 
the local standard deviation is constant throughout the 
signal, i.e., the signal becomes stationary and thus there 
is no crossover. Note the asymmetry in the value of nx 
— a much smaller value of n-^ for p — 0.05 compared 
to p = 0.95 [see Fig. ^(b-c)]. This result indicates that 
very few segments with a large standard deviation (com- 
pared to the rest of the signal) can have a strong effect 
on the anti-correlations in the signal. Surprisingly, the 
same fraction of segments with a small standard devia- 
tion (compared to the rest of the signal) does not affect 
the anti-correlations up to relatively large scales. 

For nonstationary correlated signals (a > 0.5) with 
segments characterized by two different values of the 
standard deviation, we surprisingly find no difference in 
the scaling of F{n)/n, compared to the stationary corre- 
lated signals with constant standard deviation [Fig. ^(d)]. 
Moreover, this observation remains valid for different 
sizes of the segments W and for different values of the 
fraction p of segments with larger standard deviation. 
We note that in the limiting case of very large values 
of (72 /o"! , when the values of the signal in the segments 
with standard deviation ai could be considered close to 
"zero", the results in Fig. ^(d) do not hold and we ob- 
serve a scaling behavior similar to that of the signal in 
Fig. ||(c) (see following section). 
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B. Signals with different local correlations 



Next we consider nonstationary signals which consist 
of segments with identical standard deviation (cr = 1) 
but different correlations. We obtain such signals using 
the following procedure: (1) we generate two stationary 
signals ui{i) and U2{i) (see Sec. ||) of identical length 
Nmax and with different correlations, characterized by 
scaling exponents ai and a2] (2) we divide the signals 
ui{i) and U2{i) into non-overlapping segments of size W] 
(3) we randomly replace a fraction p of the segments in 
signal ui{i) with the corresponding segments of U2{i). In 
Fig. ^(a), we show an example of such a complex non- 
stationary signal with different local correlations. In this 
Section, we study the behavior of the r.m.s. fluctuation 
function F{n)/n. We also investigate F{n)/n separately 
for each component of the nonstationary signal (which 
consists only of the segments with identical local corre- 
lations) and suggest an approach, based on the DFA re- 
sults, to recognize such complex structures in real data. 



FIG. 5. Scaling behavior of a nonstationary signal with two 
different scaling exponents, (a) Nonstationary signal (length 
Nmax = "2^^, standard deviation a = 1) which is a mixture 
of correlated segments with exponent ai =0.1 (90% of the 
signal) and segments with exponent 02 — 0.9 (10% of the 
signal). The segment size is = 20; (b) the 90% compo- 
nent containing all segments with ai =0.1 and the remain- 
ing segments (with Q2 — 0.9) are replaced by zero; (c) the 
10% component containing all segments with a2 = 0.9 and 
the remaining segments (with ai = 0.1) are replaced by zero; 
(d) DFA results for the mixed signal in (a), for the individual 
components in (b) and (c), and our prediction obtained from 
the superposition rule. 

In Fig. ^(d), we present the DFA result on such a non- 
stationary signal, composed of segments with two differ- 
ent types of local correlations characterized by exponents 
ai =0.1 and a2 = 0.9. We find that at small scales, the 
slope of F{n)/n is close to ai and at large scales, the 
slope approaches a2 with a bump in the intermediate 
scale regime. This is not surprising, since ai < a2 and 
thus F{n)/n is bound to have a small slope (ai) at small 
scales and a large slope (02) at large scales. However, 
it is surprising that although 90% of the signal consists 
of segments with scaling exponent ai, F{n)/n deviates 
at small scales (n « 10) from the behavior expected for 
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an anti-correlated signal u{i) with exponent ai [see, e.g., 
the sohd Hne in Fig. §(b)]. This suggests that the be- 
havior of F{n)/n for a nonstationary signal comprised of 
mixed segments with different correlations is dominated 
by the segments exhibiting higher positive correlations 
even in the case when their relative fraction in the signal 
is small. This observation is pertinent to real data such 
as: (i) heart rate recordings during sleep where different 
segments corresponding to different sleep stages exhibit 
different types of correlations [Q; [ii) DNA sequences in- 
cluding coding and non-coding regions characterized by 
different correlations |^,||,^ and {in) brain wave signals 
during different sleep stages [ |62[ . 

To better understand the complex behavior oi F{n)/n 
for such nonstationary signals, we study their compo- 
nents separately. Each component is composed only of 
those segments in the original signal which are char- 
acterized by identical correlations, while the segments 
with different correlations are substituted with zeros 
[see Figs. ||(b) and (c)]. Since the two components of 
the nonstationary signal in Fig. |^(a) are independent, 
based on the superposition rule [Eq. (^] , we expect that 
the r.m.s. fluctuation function F{n)/n will behave as 
^J[Fi{n)/nY + [F2{n)/nY, where Fi(n)/n and F2{n)/n 
are the r.m.s. fluctuation functions of the components in 
Fig. ||(b) and Fig. ||(c), respectively. We find a remark- 
able agreement between the superposition rule prediction 
and the result of the DFA method obtained directly from 
the mixed signal [Fig ^(d)] . This finding helps us under- 
stand the relation between the scaling behavior of the 
mixed nonstationary signal and its components. 

Information on the effect of such parameters as the 
scaling exponents ai and a2, the size of the segments W 
and their relative fraction p on the scaling behavior of the 
components provides insight into the scaling behavior of 
the original mixed signal. When the original signal comes 
from real data, its composition is a priori unknown. A 
first step is to "guess" the type of correlations (exponents 
ai and 012) present in the signal, based on the scahng be- 
havior of F{n)/n at small and large scales [Fig ||(d)]. A 
second step is to determine the parameters W and p for 
each component by matching the scaling result from the 
superposition rule with the original signal. Hence in the 
following subsections, we focus on the scaling properties 
of the components and how they change with p, a and 
W. 



1. Dependence on the fraction of segments 

First, we study how the correlation properties of the 
components depend on the fraction p of the segments 
with identical local correlations. 

For components containing segments with anti- 
correlations (0 < a < 0.5) and fixed size W [Fig. ||(b)], 
we find a crossover to random behavior (a = 0.5) at large 
scales, which becomes more pronounced (shift to smaller 



scales) when the fraction p decreases [Fig. ^(a)]. At small 
scales (n < W), the slope of F{n)/n is identical to that 
expected for a stationary signal u{i) (i.e., p = 1) with 
the same anti-correlations [solid line in Fig. ^(a)]. More- 
over, we observe a vertical shift in F{n)/n to lower values 
when the fraction p of non-zero anti-correlated segments 
decreases. We find that at small scales after rescaling 
F{n)/n by ^^p, all curves collapse on the curve for the 
stationary anti-correlated signal u{i) [Fig. ^(a)]. Since 
at small scales {n < W) the behavior of F{n)/n does 
not depend on the segment size W, this collapse indi- 
cates that the vertical shift in F{n)/n is due only to the 
fraction p. Thus, to determine the fraction p of anti- 
correlated segments in a nonstationary signal [mixture 
of anti-correlated and correlated segments. Fig. ||(a)] we 
only need to estimate at small scales the vertical shift in 
F{n)/n between the mixed signal [Fig. |5|(d)] and a sta- 
tionary signal u{i) with identical anti-correlations. This 
approach is valid for nonstationary signals where the frac- 
tion p of the anti-correlated segments is much larger than 
the fraction of the correlated segments in the mixed sig- 
nal [Fig. |5|(a)], since only under this condition the anti- 
correlated segments can dominate F{n)/n of the mixed 
signal at small scales [Fig. ||(d)]. 
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FIG. 6. Dependence of the scaling behavior of components 
on the fraction p of the segments with identical local corre- 
lations, emphasizing data collapse at small scales. The seg- 
ment size is = 20 and the length of the components is 
Nmax = 2^*^. (a) Components containing anti-correlated seg- 
ments (q = 0.1), at small scales (n < W). The slope of 
F{n)/n is identical to that expected for a stationary (p — 1) 
signal with the same anti-correlations. After rescaling F{n)/n 
by y^, at small scales all curves collapse on the curve for the 
stationary anti-correlated signal, (b) Components containing 
correlated segments (a = 0.9), at small scales (n < W). The 
slope of F{n)/n is identical to that expected for a stationary 
(p = 1) signal with the same correlations. After rescaling 
F{n)/n by at small scales all curves collapse on the curve 
for the stationary correlated signal. 



For components containing segments with positive cor- 
relations (0.5 < a < 1) and fixed size W [Fig. ||(c)], we 
observe a similar behavior for F(n)/n, with collapse at 
small scales (n < W) after rescaling by ^ [Fig. ^(b)] 
(For a > 1, there are exceptions with different rescaling 
factors, see Appendix P). At small scales the values of 
F{n)/n for components containing segments with posi- 
tive correlations are much larger compared to the values 
of F(n)/n for components containing an identical frac- 
tion p of anti-correlated segments [Fig. ||(a)]. Thus, for 
a mixed signal where the fraction of correlated segments 
is not too small (e.g., p > 0.2), the contribution at small 
scales of the anti-correlated segments to F{n)/n of the 
mixed signal [Fig. |^(d)] may not be observed, and the 
behavior (values and slope) of F{n)/n will be dominated 
by the correlated segments. In this case, we must con- 
sider the behavior of F{n)/n of the mixed signal at large 
scales only, since the contribution of the anti-correlated 
segments at large scales is negligible. Hence, we next 
study the scaling behavior of components with correlated 
segments at large scales. 

For components containing segments with positive cor- 
relations and fixed size W [Fig. ^c)], we find that at large 
scales the slope of F{n)/n is identical to that expected 
for a stationary signal u{i) (i.e., p = 1) with the same 
correlations [solid line in Fig. 0(a)]. We also observe a 
vertical shift in F{n)/n to lower values when the frac- 
tion p of non-zero correlated segments in the component 
decreases. We find that after rescaling F{n)/n by p, at 
large scales all curves collapse on the curve representing 
the stationary correlated signal u{i) [Fig. 0(a)]. Since at 
large scales (n ^ W), the effect of the zero segments 
of size W on the r.m.s. fluctuation function F(n)/n for 
components with correlated segments is negligible, even 
when the zero segments are 50% of the component [see 
Fig. 0(a)], the finding of a collapse at large scales indi- 
cates that the vertical shift in F{n)/n is only due to the 
fraction p of the correlated segments. Thus, to deter- 
mine the fraction p of correlated segments in a nonsta- 
tionary signal (which is a mixture of anti-correlated and 
correlated segments [Fig. |(a)]), we only need to estimate 
at large scales the vertical shift in F{n)/n between the 



mixed signal [Fig. ^(d)] and a stationary signal u{i) with 
identical correlations. 
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FIG. 7. Dependence of scaling behavior of components on 
the fraction p of the segments with identical local correlations, 
emphasizing data collapse at large scales. The segment size 
is = 20 and the length of the components is Nmax ~ 2'^''. 
(a) Components containing correlated segments (a = 0.9), at 
large scales (n ^ W). The slope of F{n)/n is identical to that 
expected for a stationary (p — 1) signal with the same corre- 
lations. After rescaling F{n)/n by p, at large scales all curves 
collapse on the curve for the stationary correlated signal, (b) 
Components containing anti-correlated segments (a = 0.1), 
at large scales (n 3> W). There is a crossover to random be- 
havior (a — 0.5). After rescaling F(n)/n by •y/p(l — p), all 
curves collapse at large scales. 



For components containing segments with anti- 
correlations and fixed size W [Fig. ^(b)], we find that 
at large scales in order to collapse the F{n)/n curves 
{n ^ W) [Fig. ^(a)] we need to rescale F{n)/n by 
\/p{l — p) [see Fig. 0(b)]. Note that there is a difference 
between the rescaling factors for components with anti- 
correlated and correlated segments at small [Figs, ^(a- 
b)] and large [Figs. 0(a-b)] scales. We also note that 
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for components with correlated segments (a > 0.5) and 
suffi ciently sm all p, there is a different rescaling factor 
of ^/p(^^^p) in the intermediate scale regime [see Ap- 
pendix g Fig. 

For components containing segments of white noise 
(a = 0.5), we find no change in the scaling exponent as 
a function of the fraction p of the segments, i.e., a = 0.5 
for the components at both small and large scales. How- 
ever, we observe at all scales a vertical shift in F{n)/n 
to lower values with decreasing p: F{n)/n ^ y^. 
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2. Dependence on the size of segments 



Next, we study how the scaling behavior of the com- 
ponents depends on the size of the segments W. 

First, we consider components containing segments 
with anti-correlations. For a fixed value of the fraction p 
of the segments, we study how F{n)/n changes with W. 
At small scales, we observe a behavior with a slope sim- 
ilar to the one for a stationary signal u{i) with identical 
anti-correlations [Fig. §(a)]. At large scales, we observe 
a crossover to random behavior (exponent a — 0.5) with 
an increasing crossover scale when W increases. At large 
scales, we also find a vertical shift with increasing values 
for F{n)/n when W decreases [Fig. ||(a)]. Moreover, we 
find that there is an equidistant vertical shift in F{n)/n 
when W decreases by a factor of ten, suggesting a power- 
law relation between F{n)/n and W at large scales. 
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FIG. 8. Dependence of the scaling behavior of compo- 
nents on the segment size W . The fraction p — 0.1 of the 
non-zero segments is fixed and the length of the components 
is Nmax ~ 2^". (a) Components containing anti-correlated 
segments (a — 0.1). At large scales (n ^ W), there is a 
crossover to random behavior {a — 0.5). An equidistant ver- 
tical shift in F(n)/n when W decreases by a factor of ten 
suggests a power-law relation between F{n)/n and W. (b) 
Components containing correlated segments {a = 0.9). At 
intermediate scales, F{n)/n has slope a — 0.5, indicating 
random behavior. An equidistant vertical shift in F{n)/n 
suggests a power-law relation between F{n)/n and W. 



For components containing correlated segments with a 
fixed value of the fraction p we find that in the intermedi- 
ate scale regime, the segment size W plays an important 
role in the scaling behavior of F{n)/n [Fig. ^(b)]. We 
first focus on the intermediate scale regime when both 
p = 0.1 and W — 20 are fixed [middle curve in Fig. ||(b)]. 
We find that for a small fraction p of the correlated seg- 
ments, F(n)/n has slope a — 0.5, indicating random 
behavior [Fig. |8|(b) which shrinks when p increases [see 
Appendix Fig. Thus, for components contain- 

ing correlated segments, F{n)/n approximates at large 
and small scales the behavior of a stationary signal with 
identical correlations (a — 0.9), while in the intermediate 
scale regime there is a plateau of random behavior due 
to the random "jumps" at the borders between the non- 
zero and zero segments [Fig. ||(c)]. Next, we consider the 
case when the fraction of correlated segments p is fixed 
while the segment size W changes. We find a vertical 
shift with increasing values for F{n) / n when W increases 
[Fig. 11(b)], opposite to what we observe for components 
with anti-correlated segments [Fig. H(a)]. Since the verti- 
cal shift in F{n)/n is equidistant when W increases by a 
factor of ten, our finding indicates a power-law relation- 
ship between F{n)/n and W. 
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3. Scaling Expressions 



To better understand the complexity in the scaling be- 
havior of components with correlated and anti-correlated 
segments at different scales, we employ the superposition 
rule (see and Appendix 0) . For each component we 
have 



Fin)/n = v/[-F^corr(n)/n]2 + [F,^M/n?, 



(7) 



where f con- accounts for the contribution of the 
correlated or anti-correlated non-zero segments, and 
Friindi'>T-)/n accounts for the randomness due to "jumps" 
at the borders between non-zero and zero segments in 
the component. 

Components with correlated segments (a > 0.5j 

At small scales n < W, our findings presented in Fig. ^(b) 
suggest that there is no substantial contribution from 
Fr,ndin)/n. Thus from Eq. dTl), 



F{n)/n « Fcorr(n)/" ^ boy/pn° 



(8) 



where b^n" is the r.m.s. fluctuation function for station- 
ary (p = I) correlated signals [Eq. (||) and [|6l|]. 

Similarly, at large scales n ^ W, we find that the con- 
tribution of Fi.and(«)/« is negligible [see Fig. ^(a)], thus 
from Eq. we have 



F{n)/n w i^corr /n- ^ bopn° 



(9) 



However, in the intermediate scale regime, the contribu- 
tion of F,-iin(i{n) / n to F{n)/n is substantial. To confirm 
this we use the superposition rule [Eq. (0)] and our esti- 
mates for Fcorr(n)/n at small [Eq. (||)] and large [Eq. (||)] 
scales [BSl. The result we obtain from 



^;-and(n)/n = J[F(n)/n]2 - [ho^n'-f - [b^pn^f (10) 
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FIG. 9. (a) Scaling behavior of components contain- 
ing correlated segments (a > 0.5). F(n)/n exhibits two 
crossovers and three scaling regimes at small, intermedi- 
ate and large scales. From the superposition rule [Eq. (Q)] 
we find that the small and large scale regimes are con- 
trolled by the correlations (a = 0.9) in the segments 
[Jcorr (fi) /?^ from Eqs. (0) and ^] while the intermediate 
regime [F-c^^^{n)/n ~ n from Eq. ([lo[)] is dominated by 
the random jumps at the borders between non-zero and zero 
segments, (b) The ratio F,and(Wi = 400)/i=;and(W'2 = 20) 
in the intermediate scale regime for fixed p and different val- 
ues of a, and the ratio -F'rand(Q)/frand(Qf = 0.5) for fixed p 
and W — Wi/W2- Fr^nd{n)/n is obtained from Eq. ( px| ) and 
the ratios are estimated for all scales n in the intermediate 
regime. The two curves overlap for a broad range of values 
for the exponent a, sug gest ing that -Frand (jt-) does not de- 
pend on h{a) [see Eqs. (^ij) and (p^]. 



overlaps with F{n)/n in the intermediate scale regime, 
exhibiting a slope of « 0.5: -Frand('^)/'^ ~ n'^'^ [Fig- §(&)]• 
Thus, -F'iand("-)/f^ is indeed a contribution due to the 
random jumps between the non-zero correlated segments 
and the zero segments in the component [see Fig. ^(c)]. 



Further, our results in Fig. |Kb) suggest that in the 
intermediate scale r egime F{n)/n ~ for fixed 

fraction p [see Sec. VB2|, where the power-law expo- 
nent gdct) may be a function of the scaling exponent a 
characterizing the correlations in the non-zero segments. 
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Since at intermediate scales -Frand ('T-)/'^ dominates the 
scaling [Eq. (|l0| ) and Fig. 0(a)], from Eq. (0) we find 
Fri,nd{n)/n « F{n)/n - W3^°'\ We also find that at in- 
termediate scales, F{n)Jn ^ ^/p(l — p) for fixed segment 
size W (see Appendix g Fig. |l^. Thus from Eq. (^ we 
find F^^ai{n)/n w F{n)/n ~ \/p{^ — p)- Hence we ob- 
tain the following general expression 



Fr.nd{n)/n ^ /i(a)v/p(l-p)W^^=("^n°-^ 



(11) 



Here we assume that -Frand ("■)/"■ a-lso depends directly 
on the type of correlations in the segments through some 
function h{a). 

To determine the form of gda) in Eq. (|l^), we perform 
the following steps: 

(a) We fix the values of p and a, and from Eq. ( p^ ) wc 
calculate the value of Fra,nd{n)/n for two different values 
of the segment size W, e.g., we choose Wi = 400 and 
W2 = 20. 

(b) From the expression in Eq. (p^), at the same scale n 
in the intermediate scale regime we determine the ratio: 



Fraud (W^l)/F,.and(t^2) = (W^l / )^^ • 



(12) 



(c) We plot Fi.a„d(VFi)/Fi.and(W'2) VS. a on a linear-log 
scale [Fig. |^(b)] . From the graph and Eq. ( p^ ) we obtain 
the dependence 



9c{a) 



^ log[fiand(W^l)/fiand(W^2)] 
l0g{Wl/W2) 

Ca-C/2, 0.5< a< 1 
0.50, for a > 1, 



(13) 



where C = 0.87 ± 0.06. Note that 5c(0.5) = 0. 

To determine if Fiand("-)/"- depends on h{a) in 
Eq. (pi]), we perform the following steps: 

(a) We fix the values of p and W and calculate the value 
of Frand (n) j Ti for two different values of the scaling expo- 
nent a, e.g., 0.5 and any other value of a from Eq. (10). 

(b) From the expression in Eq. (|l^), at the same scale n 
in the intermediate scale regime we determine the ratio: 



Fra„d(0.5) /i(0.5) /i(0.5) ' ^ ^ 

since gc(0.5) = from Eq. (|l3|). 

(c) We plot Fi.and(Q!)/Frand(0.5) VS. Q; On a liucar-log scale 
[Fig.ifb)] and find that when W = W1/W2 [in Eqs. (|l|) 
and (IJ)] this curve overlaps with Frand(W^i)/^rand(W2) 
vs. a [Fig. ^(b)] for all values of the scaling exponent 
0.5 < a < 1.5. From this overlap and from Eqs. ( p^ ) 
and (p^), we obtain 



h{a) 
h{0.5) 



W3' 



(15) 



for every value of a, suggesting that h(a) = const and 
thus Frand{n)/n can finally be expressed as: 



F,and(n)/n ^ Vp(1 - p)W^^='"^n°-'. (16) 



Components with anti- correlated segments (a < 0.5) 

Our results in Fig. |^(a) suggest that at small scales 
n <W there is no substantial contribution of Fi.and(")/'^ 
and that: 



F{n)/n « Fcorr("-)/"- ~ bo^n", 



(17) 



a behavior similar to the one we find for components with 
correlated segments [Eq. (||)]. 

In the intermediate and large scale regimes (n > W), 
from the plots in Fig. ^(b) and Fig. ||(a) we find the scal- 
ing behavior of F{n)/n is controlled by F^i^^din) / n and 
thus 

F{n)/n « F,a„d(n)/n ^ ^p{l - p)WS^^"'>n"-'' , (18) 

where gaia) = Ca- C/2 for < a < 0.5 [see Fig. |(b)] 
and the relation for F-caTxd{Ti) / n is obtained using the same 
procedure we followed for Eq. (pq). 



VI. CONCLUSIONS 

In this paper we studied the effects of three different 
types of nonstationarities using the DFA correlation anal- 
ysis method. Specifically, we consider sequences formed 
in three ways: {i) stitching together segments of signals 
obtained from discontinuous experimental recordings, or 
removing some noisy and unreliable segments from con- 
tinuous recordings and stitching together the remaining 
parts; {ii) adding random outliers or spikes to a signal 
with known correlations, and {in) generating a signal 
composed of segments with different properties — e.g. 
different standard deviations or different correlations. 
We compare the difference between the scaling results 
obtained for stationary correlated signals and for corre- 
lated signals with artificially imposed nonstationarities. 

(i) We find that removing segments from a signal 
and stitching together the remaining parts does not af- 
fect the scaling behavior of positively correlated signals 
(1.5 > a > 0.5), even when up to 50% of the points 
in these signals are removed. However, such a cutting 
procedure strongly affects anti-correlated signals, lead- 
ing to a crossover from an anti-correlated regime (at 
small scales) to an uncorrelated regime (at large scales). 
The crossover scale n x increases with increasing value of 
the scaling exponent a for the original stationary anti- 
correlated signal. It also depends both on the segment 
size and the fraction of the points cut out from the sig- 
nal: (1) nx decreases with increasing fraction of cutout 
segments, and (2) rtx increases with increasing segment 
size. Based on our findings, we propose an approach to 
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minimize the effect of cutting procedure on the correla- 
tion properties of a signal: When two segments which 
need to be removed are on distances shorter than the 
size of the segment, it is advantageous to cut out both 
the segments and the part of the signal between them. 

(ii) Signals with superposed random spikes. We find 
that for an anti-correlated signal with superposed spikes 
at small scales, the scaling behavior is close to that of the 
stationary anti-correlated signal without spikes. At large 
scales, there is a crossover to random behavior. For a 
correlated signal with spikes, we find a different crossover 
from uncorrelated behavior at small scales to correlated 
behavior at large scales with an exponent close to the 
exponent of the original stationary signal. We also find 
that the spikes with identical density and amplitude may 
cause strong effect on the scaling of an anti-correlated 
signal while they have a much smaller or no effect on 
the scaling of a correlated signal — when the two sig- 
nals have the same standard deviations. We investigate 
the characteristics of the scaling of the spikes only and 
find that the scaling behavior of the signal with random 
spikes is a superposition of the scaling of the signal and 
the scaling of the spikes. We analytically prove this su- 
perposition relation by introducing a superposition rule. 

(Hi) Signals composed of segments with different local 
properties. We find that 

(a) For nonstationary correlated signals comprised of 
segments which are characterized by two different values 
of the standard deviation, there is no difference in the 
scaling behavior compared to stationary correlated sig- 
nals with constant standard deviation. For nonstation- 
ary anti-correlated signals, we find a crossover at scale 
n X . For small scales n < n x , the scaling behavior is sim- 
ilar to that of the stationary anti-correlated signals with 
constant standard deviation. For large scales n > rix, 
there is a transition to random behavior. We also find 
that very few segments with large standard deviation can 
strongly affect the anti-correlations in the signal. In con- 
trast, the same fraction of segments with standard devi- 
ation smaller than the standard deviation of the rest of 
the anti-correlated signal has much weaker effect on the 
scaling behavior — rix is shifted to larger scales. 

(b) For nonstationary signals consisting of segments 
with different correlations, the scaling behavior is a su- 
perposition of the scaling of the different components 
— where each component contains only the segments 
exhibiting identical correlations and the remaining seg- 
ments are replaced by zero. Based on our findings, we 
propose an approach to identify the composition of such 
complex signals: A first step is to "guess" the type of cor- 
relations from the DFA results at small and large scales. 
A second step is to determine the parameters defining the 
components, such as the segment size and the fraction of 
non-zero segments. We studied how the scaling charac- 
teristics of the components depend on these parameters 
and provide analytic scaling expressions. 
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APPENDIX A: SUPERPOSITION RULE 

Here we show how the DFA results for any two signals 
/ and g [denoted as Ff{n) and Fg{n)] relate with the DFA 
result for the sum of these two signals f + g [denoted as 
Ff+g{n), where n is the box length (scale of analysis)]. 
In the general cases, we find \Ff — Fg\ < Ff+g < Ff + Fg. 
When the two signals are not correlated, we find that the 



following superposition rule is valid: Fj^^ ~ f 



F^ + Fl 



Here we derive these relations. 

First we summarize again the procedure of the DFA 
method [^. It includes the following steps: starting with 
an original signal u{i) of length N^axi we integrate and 

fe 

obtain y{k) = ^ ("(i) ~ where {u) is the mean of 

u{i). Next, we divide y{k) into non-overlapping boxes of 
equal length n. In each box we fit the signal y{k) using a 
polynomial function ?/„(fc) — ao-l-aix(fc)-|-a2a;^(fc)-|-...-|- 
asx'^ik), where x{k) is the x coordinate corresponding to 
the kth signal point. We calculate the r.m.s. fluctuation 



function F{n) = 



N„ 



N„ 



E [y{k)-y,,{k)Y. 



To prove the superposition rule, we first focus on one 
particular box along the signal. In order to find the an- 
alytic expression of best fit in this box, we write 

n 

/(ao,...,a«) =^[y(fc)- (ao + ...-l-a,a;'^(A:))]2, (Al) 
fe=i 

where am, m = 0, s are the same for all points in this 
box. "Best fit" requires that a™, m — 0, s satisfy 



dl 

dam. 



0, TO = 0, ...s 



(A2) 



Combining Eq. (Al) with (A2) wc obtain s + 1 equations 
ym = aotmo + aitmi--- + agtms, TO = 0, s (A3) 
where 

n n 

ym^Y. V{k)x^\k),tm, = J2 = 0, s (A4) 



fc=i 



k=l 



From Eqs. (A3) we determine ag, oi, , 



14 



For the signals /, g and f + g after the integration, in 
each box we have 

fm = flO^mO + ai^ml--- + C-stms, TO = 0, S 
g-m = Oo^mO + a'l^ml--- + a'sUns, m — 0, S 

if + 9)m = flo^mO + a'itrni--- + a'^Uns, m = 0, s (A5) 



where fm, gm and {f+g)m correspond to inEqs. (AS). 

Comparing the three groups of equations in Eqs. (^), 
we find that, when we add the first two groups together, 
the left side becomes /m + 3m = (/ + .9)m , which is pre- 
cisely the left side of the third group of equations. Thus 
we find 



flm + «mj m = 0, . 



(A6) 



Yg{k) — where Yf{k) and Yg{k) fluctuate around 



fc=i 



zero. When Yf{k) and Fgffc) are not correlated, the value 
of the third term in Eq. ( [A9| ) is close to zero and we ob- 
tain the following superposition rule 



F 



f+g 



Fi + f: 



(A12) 



APPENDIX B: STRONGLY CORRELATED 
SEGMENTS 



and for each point k in every box, the polynomial fits for 
the signals /, g and f + g satisfy 



(/ + 3)n(fc) =/„(fc)+5n(fc). 



(A7) 



This result can be extended to all boxes in the signals. 
For the signal f + g we obtain 



1 



N„ 



N„ 



E ifik)-fn{k)]' + [9ik)-gnik)r 



k=l 



-2[/(fc)-/„(fc)][g(fc)-5„(fc)]. 



(A8) 



After the substitutions f{k) — fn{k) = Yf{k) and g{k) 
gn(k) — Ygik), we rewrite the above equation as 



k=l 



k=l 



N„ 



+2 J2 Yf{k)Yg(k) 
2 



k=l 



N„ 



Fi + Ft 



N„ 



(A9) 



k=l 



In the general case, we can utilize the Cauchy inequal- 



ity 



N„ 



E ^/wn(fc) 



fc=i 



1/2 



1/2 



< E (^/W)' E (^.w)' 



(AlO) 



\ k=l 

and we find 



fc=i 



{Ff - Fgf < Ff^g < {Ff + Fgf 
\Ff - Fg\ < Ff+g < Ff + Fg. 



(All) 



From Eqs. ( [A3| ) for m = 0, in every box we have 
y{k) = Vn{k). Thus we obtain Yf{k) = 

k=l k=l k=l 



For components containing segments with stronger 
positive correlations (a > 1) and fixed W — 20, we find 
that at small scales {n < W), the slope of F{n)/n does 
not depend on the fraction p and is close to that ex- 
pected for a stationary signal u{i) with identical corre- 
lations (Fig. [To| ). Surprisingly we find that in order to 
collapse the F{n)/n curves obtained for different values 
of the fraction p, we need to rescale F{n)/n by y^pil — p) 
instead of which is the rescaling factor at small 
scales for components containing segments with corre- 
lations a < 1. Thus a = 1 is a threshold indicating 
when the rescaling factor changes. Our simulations show 
that this threshold increases when the segment size W 
increases. 

For components containing a sufficiently small fraction 
p of correlated segments (a > 0.5), we find that in the 
intermediate scale regime there is a crossover to random 
behavior with slope 0.5. The F{n)/n curves obtained 
for different values of p collapse in the intermediate scale 
regime if we rescale F{n)/n by y^p{l — p) (Fig. |l^). We 
note that this random behavior regime at intermediate 
scales shrinks with increasing the fraction p of correlated 
segments and, as expected, for p close to 1 this regime 
disappears (see the p = 0.9 curve in Fig. p^). 
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FIG. 10. Dependence of the scaling behavior of compo- 
nents on the fraction p of the segments with strong positive 

correlations {a = 1.2). The segment size is = 20 and 
the length of the compon ents is Nmax ~ 2^". After rescal- 
ing F{n)jn by 1^/^(1 — p), all curves collapse at small scales 
(n < W) with slope 1.2 and at intermediate scales with slope 
0.5. The intermediate scale regime shrinks when p increases. 
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